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IV 


ABSTRACT 


Silicon  Blocked  Impurity  Band  (BIB)  detectors  are  state-of-the-art  devices 
to  detect  light  in  the  near  to  mid  infrared  range  (5-40  //m).  Numerical  modeling 
of  BIB  detectors  is  performed  using  a  four-region  finite  difference  approach  to 
study  the  role  of  space  charge  in  C-V  (capacitance-voltage)  profiling  of  minority 
carrier  doping  and  the  role  of  blocking  layer  thickness  and  minority  doping 
concentration  in  an  alternative  bias  scheme.  Compensation  in  the  blocking  layer 
is  found  to  play  a  critical  role  in  determining  the  net  voltage  drop  in  this  part  of  the 
device  under  alternate  polarity  bias.  The  effect  of  space  charge  at  the  blocking 
layer/active  layer  interface  on  the  measured  low  temperature  C-V  distribution  is 
modeled  as  a  function  of  the  doping  interface  between  the  two  layers.  The 
magnitude  of  the  space  charge  can  cause  large  deviations  in  the  measurement 
of  minority  doping  concentration  from  the  idealized  case  which  assumes  a  space- 
charge  free  blocking  layer  and  interface.  Accurately  determining  these  minority 
doping  concentrations  is  a  crucial  step  toward  solving  material  growth  challenges 
in  proposed  far  infrared  Ge  and  GaAs  devices. 
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I.  INTRODUCTION 


A.  BLOCKED  IMPURITY  BAND  (BIB)  DETECTORS 

Blocked  Impurity  Band  (BIB)  detectors,  also  known  as  Impurity  Band 
Conduction  (IBC)  detectors,  are  highly  efficient,  radiation  hardened  extrinsic 
photodetectors  currently  used  to  detect  light  in  the  5-40  //  m  range.  They  were 
first  developed  in  1986  by  Petroff  and  Stapelbroek  at  the  Rockwell  International 
Science  Center  [1], 

BIB  detectors  consist  of  two  layers  fabricated  through  liquid  or  vapor 
phase  epitaxy.  One  region,  the  active  layer,  is  heavily  doped  and  is  responsible 
for  photon  absorption  through  donor  (or  acceptor)  ionization.  This  high  doping 
level  allows  for  greater  efficiency  and  reduced  volume  and  is  responsible  for  BIB 
replacement  of  standard  photoconductors  in  many  applications.  The  second 
region  is  a  highly  pure  layer  called  the  blocking  layer.  This  region  limits  the 
current  to  carriers  excited  to  the  conduction  (or  valence)  band  of  the  device. 

The  small  energies  used  for  detection  vary  with  the  material  and  dopant 
used.  For  an  As  doped  Si  device,  the  ionization  energy  is  54  meV.  This 
translates  into  a  cutoff  wavelength  of  about  26  //m  and  requires  an  operating 
temperature  of  4-10  K.  This  wavelength  response  has  been  extended  with  Si:Sb 
BIBs  which  have  ionization  energies  of  43  meV  and  wavelength  extension  out  to 
37  //m  [2], 

Current  work  is  focusing  on  extending  BIB  response  into  the  far-infrared 
range  using  Ge:Ga  and  GaAs.  Ge:Ga  has  ionization  energies  of  10  -  11  meV 
with  wavelength  cutoff  response  to  100  //  m,  while  uniaxial  compression  can 
lower  these  ionization  energies  and  extend  response  out  to  200  //m.  GaAs 
doped  with  Si  or  Te  demonstrates  ground  state  ionization  energies  of  5.7  meV 
with  the  first  excited  state  transition  at  an  energy  of  4.3  meV.  These  energies 
represent  wavelength  response  at  217  and  286  //  m  respectively. 


1 


Due  to  the  small  energies  involved  in  detecting  light  in  the  infrared  and 
especially  far-infrared  range,  the  device  must  be  cooled  to  prevent 
donor/acceptor  sites  from  becoming  ionized  producing  excessive  thermal- 
excitation  into  the  conduction/valence  band.  Examples  of  the  various  materials 
and  wavelength  responses  described  can  be  seen  in  Figure  1.  The  longer 
wavelength  response  materials  offer  great  potential,  however  fabrication  issues 
continue  to  hinder  the  realization  of  these  devices  [3], 


Figure  1  BIB  wavelength  extension. 


Foremost  among  these  challenges  is  growing  the  extremely  thin,  highly 
pure  blocking  layer  the  devices  require.  It  can  be  difficult  to  grow  blocking  layers 
to  the  purity  levels  desired.  First,  the  majority  doping  level  ideally  must  be  very 
low  (-1011  cm'3).  In  addition,  the  minority  doping  must  be  even  lower  relative  to 
the  majority  doping  level.  The  compensation  of  a  device  region,  defined  as 
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must  be  controlled  precisely  and  kept  very  low  for  many  applications.  Even  when 
these  purity  levels  are  attainable,  the  layer  also  needs  to  be  grown  very  thin,  and 
often  the  blocking  layer  never  attains  the  ideal  doping  level.  In  these  cases  the 
blocking  layer  is  simply  a  transition  area  from  the  heavily  doped  active  layer. 
Diffusion  of  dopants  from  the  active  layer  into  the  blocking  layer  adds  to  this 
graded  region.  The  end  result  of  these  combined  factors  is  a  transition  area  with 
enough  doping  to  allow  significant  dark  current. 

The  difficulty  of  these  growth  challenges  leads  naturally  to  application  of 
various  methods  of  device  characterization.  Inability  to  directly  measure  doping 
characteristics  only  makes  growth  challenges  more  difficult  to  solve,  and 
investigation  and  analysis  of  one  of  these  key  techniques,  low  temperature  C-V 
profiling  for  minority  doping  concentrations,  will  be  analyzed  in  detail. 

B.  MILITARY  RELEVANCE 

BIBs  were  originally  intended  to  provide  epitaxial  style  detectors  for  the 

mid-IR  range  to  replace  conventional  bulk  Si  photoconductors,  which  were  highly 

sensitive  to  degenerative  radiation  effects.  The  primary  motivation  was  for 

defense  applications,  specifically  for  ballistic  missile  defense.  The  IR  response 

of  these  detectors  makes  them  ideal  for  cold  body  detection  and  tracking,  an 

important  feature  for  re-entry  and  late-stage  interception  and  destruction  of  a 

missile  [4],  In  addition  a  great  deal  of  spectroscopy  information  is  available  in 

this  range,  so  chemical  identification  using  these  devices  is  of  significant  interest. 

The  highly  touted  tera-hertz  frequencies  also  lie  within  this  range,  which  has 

potential  for  significantly  outperforming  x-ray  and  metal  detector  devices  in 

security  fields.  In  this  area,  infrared  spectroscopy  is  of  interest  to  identify 

chemical  species  in  closed  containers.  However,  these  highly  sensitive  detector 

arrays  have  always  found  a  niche  in  space-based  mid-IR  astronomy.  Used  for 

deep  space  imaging  and  spectroscopy,  Si:As  and  Si:Sb  BIB  devices  are  currently 

on  the  Spitzer  Space  Telescope  and  a  Si:As  BIB  array  is  planned  for  the  James 

Webb  Space  Telescope  to  be  launched  in  201 1  [2,5], 
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C.  BIB  OPERATION 

A  complete  description  of  BIB  operation  is  difficult  because  many  different 
doping  schemes  and  bias  configurations  exist.  Several  principles  are,  however, 
common  to  all  BIB  detectors.  The  operation  of  the  BIB,  regardless  of  type,  can 
be  tied  to  two  key  characteristics:  the  use  of  a  heavily  doped  active  layer  for 
collection  of  light  and  a  pure  blocking  layer  which  prevents  flow  of  the  hopping 
current  along  the  dopant  band.  All  BIBs  operate  by  establishing  a  depletion 
region  in  the  active  layer  which  collects  light  through  ionization  of  dopant  sites. 
Once  excited  to  the  conduction  band  (or  in  the  case  of  holes  to  the  valence 
band),  these  carriers  create  a  photocurrent. 

In  standard  bias  configuration,  a  bias  is  applied  to  the  blocking  layer 
contact,  thus  depleting  through  the  highly  pure  blocking  layer  into  the  active  layer 
where  the  necessary  depletion  region  is  established.  The  simplest  starting  point 
is  shown  in  Figure  2  with  an  n-type  BIB  such  as  Silicon  doped  Antimony  (Si:Sb), 
in  which  the  majority  dopant  has  an  extra  electron,  and  thus  constitutes  a  donor 
site.  All  of  the  minority  acceptor  sites  in  the  active  layer  are  ionized  by  the 
majority  dopants.  In  this  schematic  the  electric  field  in  the  device  is  represented 
with  a  red  line.  A  positive  bias  is  applied  to  the  blocking  layer  contact,  which 
depletes  into  the  active  layer. 

Due  to  the  low  temperature  of  operation  the  electrons  remain  thermally 
within  the  donor  band,  but  are  free  to  move  along  it  in  the  active  layer.  Electrons 
will  move  to  occupy  the  ionized  donors,  leaving  a  net  space  charge  of  ionized 
minority  dopants  ( N~ ).  The  depletion  width  is  therefore  determined  by  the 

concentration  of  minority  dopants.  In  the  case  of  an  n-type  device  this  means  a 
positive  bias  must  be  applied  in  order  to  terminate  the  electric  field  lines  at  the 
negative  ionized  acceptors  {NA).  When  light  ionizes  a  donor  site  in  the  depletion 
width,  the  bias  moves  free  electrons  toward  the  blocking  layer  contact  and  drives 
ionized  donor  (N^)  sites,  the  hopping  current,  toward  the  active  layer  contact. 
Thus,  the  photocurrent  is  collected  in  a  standard  biased  n-type  BIB. 


4 


Figure  2  Schematic  diagram  of  a  standard  biased  n-type  BIB.  Red  line 

indicates  the  electric  field  distribution. 

P-type  BIBs  such  as  galluim  doped  germanium  (Ge:Ga)  in  standard  bias 
configuration,  Figure  3,  operate  similarly  only  the  majority  dopants  are  acceptors 
rather  than  donors.  Once  again,  the  minority  dopants,  in  this  case  ND ,  are 

compensated  by  the  majority  dopants.  The  ionized  donor  sites  {N+D)  create  a 

positive  depletion  width  which  requires  applying  a  negative  bias  to  the  blocking 
layer  contact.  When  the  acceptors  are  ionized  by  incident  light,  the  resulting  hole 
is  pulled  along  the  valence  band  to  the  blocking  layer  contact,  and  the  NA  carrier 

is  pushed  to  the  active  layer  contact.  Again,  this  is  the  photocurrent  collection 
method  in  a  standard  biased  p-type  BIB  [6,7], 
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Figure  3  Schematic  diagram  of  standard  biased  p-type  BIB.  Red  line 

indicates  the  electric  field  distribution. 

These  two  relatively  simple,  oppositely  doped  devices  can  result  in  a  great 
deal  of  complexity,  especially  with  the  variety  of  currents  and  space  charge  that 
must  be  included  in  any  analysis.  The  photogenerated  majority  carriers, 
dependent  on  the  majority  dopant,  create  the  most  obvious  source  of  current, 
however  the  ionized  majority  dopant  sites  also  create  the  hopping  current.  The 
minority  dopant  brings  another  charge  to  be  considered,  specifically  for  depletion 
regions,  and  drift  and  diffusion  currents  between  the  two  layers  and  from  the 
contacts  must  be  included.  The  effects  of  these  various  charges  lead  to 
complex,  sometimes  counter  intuitive  results.  Analysis  of  space  charge  effects, 
particularly  in  the  blocking  layer,  on  characterization  techniques  will  be  analyzed 
in  depth  in  this  thesis. 
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II.  EXPERIMENTAL  SETUP 


A.  COMPUTATION  SETUP 

All  of  the  research  into  BIBs  at  the  Naval  Postgraduate  School  is 
computational  modeling  and  device  simulation.  The  computational  model 
outputs  13  different  variables  for  each  iteration,  and  a  single  run  generally 
requires  fifteen  to  twenty  thousand  iterations  depending  upon  the  applied  bias. 
This  research  required  over  500  simulations  which  would  not  have  been  possible 
without  several  computationally  powerful  computers.  The  BIB  lab  consists  of 
three  machines  with  AMD  (3000+)  Athlon  64  Bit  processors,  Suse  Linux  9.1 
operating  systems,  and  2  GB  400  DDR  RAM.  Two  of  the  computers  have  80  GB 
hard  drives,  while  the  third  has  a  120  GB  hard  drive.  Finally,  the  data  analysis 
was  performed  on  a  fourth,  Windows  XP  machine  using  Sigma  Plot. 

B.  FINITE  DIFFERENCE  COMPUTER  MODEL 

The  three  Linux  computers  run  a  finite  difference,  one  dimensional  model 
which  calculates  steady  state  values  for  the  electric  field,  carrier  concentrations, 
free  carrier  current,  and  ionized  impurity  band  (hopping)  current  as  functions  of 
position.  This  model  has  been  used  previously  to  study  field  distributions  in  BIB 
devices  [8,9],  The  model  arrives  at  these  parameters  by  solving  Poisson’s 
equation,  current  continuity,  and  detailed  balance  in  one  dimension.  This  model 
is  written  to  solve  these  equations  for  a  p-type  device,  although  it  is  simple  to 
model  n-type  devices  through  sign  adjustments  to  the  model  outputs. 

Using  Poisson’s  equation 

££o^  =  q(p-n  +  Nd~Na)  (2-1 ) 

it  is  possible  to  solve  for  the  electric  field  E  as  a  function  of  space  charge.  Here 
ss0  is  the  permittivity  of  the  region,  p  and  n  are  the  concentrations  of  free 
electrons  and  holes  at  each  position,  Nd+  is  the  concentration  of  ionized  donor 
sites,  and  NA  is  the  concentration  of  ionized  acceptor  sites.  The  ionized  minority 
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population  (for  a  p-type  device  )  is  treated  as  the  minority  doping  ND  in  the 

model  due  to  compensation  by  the  majority  dopant.  In  the  model,  two  mobile 
charge  carriers  are  included  in  calculations,  the  majority  carrier  (for  a  p-type 
device  free  holes  in  the  valence  band)  and  the  ionized  majority  dopant  (ionized 
sites  in  the  dopant/acceptor  band). 

The  current  continuity  equations  for  positive  and  negative  ionized  charge 
carriers  are 


dp 

dt 


G-R 


(2.2) 


dn 

dt 


=  R-G 


(2.3) 


In  these  equations,  G  represents  the  generation  rate  and  R  represents  the 
recombination  rate  of  the  free  carriers.  The  generation  rate  can  be  defined  as 

G  =  (r  +  ovpPl)(NA-N-)  (2.4) 


in  which  y  is  the  optical  generation  of  free  carriers  and  the  second  term 
represents  the  thermal  generation  rate.  This  second  term  comes  from  a  detailed 
balance  argument  which  requires  the  thermal  emission  rate  from  remaining 
neutral  acceptor  sites  to  equal  the  recombination  current. 

R  ~  avpPN A  (2.5) 

In  relating  these  two  equations,  a  is  the  capture  cross  section  of  an  ionized 
acceptor  site,  v  is  the  effective  speed  of  the  ionized  majority  free  carrier,  and  pi 

is  an  unknown.  To  solve  for  this  value,  a  new  factor  a  must  be  solved  for  in  the 
relation  between  the  thermal  emission  rate  and  the  recombination  current. 

«(na-na)  =  ^pPna  (2-6) 


Here,  at  thermal  equilibrium 


-qEF 

p  =  Nre  kT 
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(2.7) 


in  which  Nv  is  the  density  of  states  in  the  valence  band  and  Ef  is  the  energy  of 
the  Fermi  level  above  the  valence  band.  The  thermal  equilibrium  relationship 
between  ionized  and  neutral  acceptors  is  the  Fermi  relationship 


This  allows  one  to  relate  the  Generation  and  Recombination  rates 

G-R  =  (r  +  (TupP] )(na  -N~)~  aupPN~  (2.11) 

which  must  be  zero  for  y  equal  to  zero  per  the  assumptions  made  for  the 
calculation. 

Referring  back  to  the  continuity  equations  (Equations  2.2  and  2.3) 
equations  for  Jn  and  Jp  are  also  required.  These  currents  are  the  combination  of 
both  diffusion  and  drift  currents  for  electrons  and  holes. 


dn 

Jn  =  qD, ,  IT  +  WMx)E (*) 
ox 

(2.12) 

dp 

JP  =(lDp-^  +  ?MpP(x)e(x) 

(2.13) 

The  drift  components  of  these  currents  include  the  hopping  current  of  the  region, 
and  to  solve  for  this  the  model  uses  an  analytical  expression  for  the  mobility  of 
these  ionized  acceptors.  The  model  assumes  the  simplified  case  of  uniformly 
spaced  dopants 
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(2.14) 


3.5  —x 

/u  =  jU0x  e 

x  =  2  [Nfa,']'  (2.15) 

with  mobility  treated  as  a  function  exponentially  dependent  on  the  Bohr  radius  ae, 
the  majority  dopant  concentration  N A ,  and  a  temperature  dependent  pre-factor 

A,  [10]. 


Position  (micron) 

Figure  4  Plot  of  doping  profile  as  a  function  of  position  for  four  values  of 

grade  parameter  (g  from  equation  2.16). 

Finally,  in  defining  the  entire  device,  the  model  includes  smooth  transitions 

between  layer  dopings  through  the  use  of  an  empirical  hyperbolic  function. 


A^  +  OV-AO 

i +ew-m 


(2.16) 


Here,  the  doping  concentration  at  any  given  point  on  an  interface  is  a  function  of 
the  doping  Ni  in  the  previous  region,  the  doping  level  N2  in  the  next  region,  the 
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position  a  of  the  interface,  the  current  position  x,  and  the  grading  parameter  g 
which  is  set  prior  to  running  the  program.  Figure  4  shows  the  effects  varying  this 
grade  parameter  has  on  the  doping  concentration  transition  between  different 
layers  [8], 

Employing  this  powerful  set  of  equations  requires  an  extensive  input  of 
variables  to  define  the  characteristics  of  the  BIB  device  to  be  modeled.  Figure  5 
shows  an  example  of  these  input  files  with  many  of  the  key  parameters  labeled. 


Run  Number 


Doping  Grade  (cm) 


End  Boundary: 
Active  Layer  (cm) 


End  Boundary: 
Blocking  Layer 
(cm) 


l.eigth 
0(  25 


1100 
Edges 
3 

Edgeposit| 
0.000025 
0.002  — I 
0.002475 
Mu  Mue(Pre-f) 
3E6  6 


Windotw 
1  8e-4 
Grade 
0.00001 


Blocking  Layer 
Acceptors  (cm  3) 


Cell  max 
1.0E-6 
nl 

1.00E+17 


Cellmiji 

3.0E-8 

n2 

1.00E-16 


0.000002—1  1.00E+17 
0.00001  1.00E+17 


BOHR 

15E-8 


GAMMA 

5.0e-8 


Active  Layer 
Acceptors  (cm*3) 

Blocking  Layer 
Donors  (cm  3) 

Active  Layer 
Donors  (cm-3) 

Bias  DdtaposC  Tim  sup 
2.0000  ).0  33' 

NA  ND 

2.00E+17  1 .00E+13 

— 1.00E+11 
1.00E+12 

DGAMMA  Ea  SVP 
0.0  0.043  2E-5 


1  00E-16  I — 1.00E+12 
1.00E-16  1.00E+13 


Temperature  (K) 


Bohr  Radius  (cm) 

Active  Layer  Optical 

Hopping  Mobility  Generation  (s1) 

, ,  (cm-  V  s) 

Active  Layer 
Mobility  Ictn-  V  s) 


Ionization 
Energy  (eV) 

Sigma  Velocity 
Product  (cm3/s) 


Figure  5  BIB  input  file 

The  model  treats  the  BIB  as  a  four  layer  device,  so  all  of  these  layers  must 
be  defined  in  the  input  file.  The  two  contacts  are  treated  as  heavily  doped  ohmic 
regions  which  supply  free  carriers  at  all  temperatures.  In  a  p-type  device  the 
model  treats  the  BIB  as  a  system  with  free  holes  and  ionized  acceptors,  so  the 
contact  is  assumed  to  supply  holes.  Usually  the  contacts  are  assigned 
thicknesses  of  0.25  //m.  As  seen  in  Figure  5  the  active  layer  and  the  blocking 
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layer  can  also  assume  whatever  thickness  is  desired.  The  majority  and  minority 
dopant  concentrations  are  assigned  to  each  region,  and  the  grades  between 
these  values  can  be  assigned  for  all  three  interfaces  in  the  device. 

The  input  file  also  allows  for  defining  different  host  materials,  and  these 
variables  include  the  mobility  of  free  carriers,  the  mobility  of  ionized  acceptors 
(ionized  majority  dopants),  the  Bohr  radius  aB,  the  majority  dopant  ionization 
energy,  and  the  cross  section  velocity  product  crt>p.  All  of  these  material 

parameters  are  listed  on  the  bottom  line  of  the  input  data  file. 

Finally,  the  input  file  includes  individual  run  specific  variables  such  as  the 
operation  temperature  in  Kelvin,  the  Optical  Generation  rate  y  which  is  a 
measure  of  the  flux  on  the  device,  and  the  bias  applied  to  the  device. 

The  model  takes  this  input  file  of  the  BIB  characteristics  and  calculates  the 
13  variables  shown  in  Figure  6.  It  arrives  at  this  solution  by  stepping  up  the 
applied  bias  from  zero  to  the  input  file  value  in  steps  of  0.01  mV.  It  arrives  at  a 
steady  state  solution  for  each  position  at  a  set  bias,  and  then  using  the  previous 
solution  and  new  values  solves  for  the  next  slightly  higher  bias.  The  final  output 
file  is  the  value  of  all  the  variables  at  each  incremental  position  in  the  device,  for 
the  final  bias  provided  on  the  input  file. 


Figure  6  BIB  output  file 
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Achieving  a  timely  solution  in  BIB  modeling  is  often  dependant  upon  how 
the  problem  is  defined.  Each  increment  of  0.01  mV  requires  a  large  number  of 
iterations  as  the  program  steps  through  the  device.  The  length  of  the  device  and 
complexity  (doping  concentrations  and  grade  parameters)  will  determine  how 
large  or  small  the  total  number  of  iterations  will  be.  Likewise,  the  lower  the 
temperature  and  the  optical  flux  values,  the  more  iterations  required.  Simulations 
for  Si  BIBs  tend  to  run  significantly  faster  than  those  for  GaAs  or  Ge  BIBs  due  to 
more  ideal  field  distributions.  The  end  result  is  that  a  solution  can  be  completed 
in  a  matter  of  minutes  or  in  weeks.  Usually  approximations  can  be  made  to 
reduce  the  time  required  to  complete  a  simulation,  however  the  model  allows  for 
highly  accurate,  and  time  consuming,  calculations  when  the  need  exists. 


13 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 


14 


III.  C-V  PROFILING  METHOD 


A.  GROWTH  CONSTRAINTS 

Challenges  of  growing  extremely  thin,  high  purity,  and  low  compensation 
materials  are  currently  limiting  the  extension  of  BIBs  into  the  far-infrared  regime. 
Minority  doping  levels  in  the  active  layer  are  of  critical  importance  and  pose  a 
special  challenge  in  the  growth  of  compound  semiconductors  such  as  GaAs. 
Because  the  depletion  width  in  the  active  layer  is  a  function  of  the  compensation, 
the  depletion  widths  become  prohibitively  small  as  compensation  levels  get  too 
high.  This  can  be  seen  in  Figure  7,  which  compares  two  GaAs  BIBs,  one  with 
0.1%  compensation  and  the  other  with  10%  compensation  in  the  active  layer. 


Figure  7  Active  layer  compensation  effects  in  GaAs  (solid  line  0.1%, 

dashed  line  10%  compensation) 
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For  the  case  with  a  compensation  of  0.1%,  the  device  would  function  as  a 
BIB  photodetector  because  the  depletion  width  extends  into  the  active  region  due 
to  the  lower  compensation.  As  a  result  the  photocurrent  from  ionized  carriers  in 
the  depletion  region  could  be  collected.  The  second  compensation  of  10%  is  so 
large  the  depletion  region  shrinks  to  nothing,  and  any  photo-ionized  carriers  in 
the  active  layer  will  recombine  rather  than  be  collected  as  photocurrent. 

Difficulties  in  achieving  extremely  thin,  pure  blocking  layers  also  raise 
minority  doping  issues.  The  blocking  layer  ideally  has  an  extremely  low  minority 
doping  level  to  minimize  space  charge  effects.  In  addition  compensation  should 
be  minimized  such  that  as  little  field  is  dropped  in  the  layer  as  possible.  This 
material  constraint  can  be  difficult  to  attain,  and  in  the  case  of  Ge  significant 
grading  in  the  blocking  layer  active  layer  interface  requires  a  thicker  layer  to  allow 
the  doping  levels  to  transition  to  the  required  purity  level.  Figure  8  illustrates 
these  difficulties  as  the  depletion  width  decreases  significantly  with  both  thicker 
blocking  layers  and  with  higher  compensation  levels. 


16 


Figure  8  Electric  field  profiles  for  varied  blocking  layer  thicknesses  (10  and 

20  /j  m)  and  minority  doping  levels  in  a  Ge  BIB. 

Common  to  both  of  these  situations  is  the  importance  of  knowing  and  then 
controlling  compensation  levels  throughout  the  BIB  device  for  the  development  of 
longer  wavelength  response  BIBs.  Certainly  the  difficulty  of  controlling 
compensation  levels  can  be  attributed  to  many  factors,  but  the  assumptions  and 
error  in  measurement  techniques  adds  to  an  already  complex  problem  [1 1,12]. 

B.  C-V  PROFILING  METHOD  THEORY 

In  isolated  bulk  samples  or  epitaxial  layers,  compensation  levels  are 
generally  measured  using  variable  temperature  Hall  effect.  The  minority  doping 
concentration  affects  the  statistics  of  recombination  and  causes  variations  in 
slope  of  the  majority  carrier  “freeze-out”  curve.  However,  this  technique  can 
generally  not  be  applied  to  layers  within  a  multi-layer  structure.  A  different 
approach  is  therefore  required  to  measure  BIB  material  parameters  within  the 
device  structure. 
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The  current  technique  used  to  measure  compensation  in  Blocked  Impurity 
Band  Detectors  is  a  low  temperature  application  of  the  C-V  profiling  method.  C-V 
profiling  involves  measuring  the  change  in  capacitance  of  the  active  layer  in  a 
BIB  detector  for  a  very  small  change  in  applied  bias.  An  initial  bias  is  applied  to 
the  device  such  that  the  depletion  width  extends  into  a  region  where  the  doping 
characterization  is  desired.  A  slight  voltage  variation  is  applied,  and  the  change 
in  capacitance  of  the  device  is  measured.  Using  this  AC  and  the  small  change 
in  voltage,  the  doping-related  space  charge  in  the  region  can  be  calculated  [13]. 

In  standard  C-V  profiling,  we  assume  that  the  minority  dopant  sites  are 
compensated  by  the  majority  dopants  and  the  temperature  of  the  device  is  high 
enough  to  allow  free  movement  of  the  majority  carriers.  When  the  depletion 
width  is  established  for  this  case,  the  total  space  charge  remaining  is  a 
combination  of  both  the  ionized  majority  and  minority  dopant  sites,  Nmaj-Nmm. 

The  majority  concentration  is  significantly  higher  than  the  minority  concentration, 
so  this  method  is  used  to  measure  net  majority  doping  in  devices.  For  the  low 
temperature  variation  the  minority  dopants  are  once  again  compensated,  but  the 
low  thermal  energy  causes  the  majority  carriers  to  be  bound  to  their  dopant  sites. 
The  result  is  that  the  space  charge  in  the  depletion  region  arises  from  only 
minority  dopants. 

This  low  temperature  C-V  profiling  method  begins  with  the  approximation 
that  within  the  depletion  width  the  space  charge  is  proportional  to  the  minority 
doping  concentration  Nd 

P  =  qND  (3.1) 


From  this  approximation  the  depletion  width  can  be  calculated. 


W  = 


Vbi 


v 


V 


kT '' 
C1  , 


(3.2) 


In  this  equation  Vbi  is  the  built-in  potential,  V  is  the  applied  potential,  and  the 

kT 

term  —  takes  into  account  the  majority-carrier  distribution  tail.  Multiplying  the 

q 
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space  charge  with  the  depletion  width,  Equation  3.1  with  3.2  gives  the  total  space 
charge  in  the  depletion  region. 


Qsc  =  J2£sqNL 


(  IT  ^ 


V 


<1 ) 


(3.3) 


The  derivative  of  Equation  3.3  with  respect  to  change  in  bias  voltage  allows  for 
the  calculation  of  the  capacitance. 


C  = 


<i£snd 


v  -V - 

y  bi  v 


kT 


(3.4) 


The  desired  parameter  is  the  minority  concentration,  so  Equation  3.4  can  be 
manipulated 


d(\/C2)_  2 

dV  qssND 


(3.5) 


such  that  the  minority  concentration  can  be  isolated. 


Nr 


2 


1 

d(l/C2)/dV 


(3.6) 


Experimental  determination  of  ND  requires  two  measurements  of  capacitance  at 

slightly  different  applied  biases.  This  surprisingly  simple  solution  hides  several 
important  assumptions.  The  starting  point  of  the  calculations  is  the  assumption 
that  all  the  space  charge  is  in  the  active  layer,  and  that  the  space  charge  consists 
exclusively  of  ionized  minority  sites.  Most  significant  of  these  assumptions  is  that 
the  blocking  layer  purity  is  high  enough  to  ignore  any  space  charge  in  the  layer. 

The  goal  of  this  work  is  to  perform  numerical  simulations  of  the  C-V 
method  that  include  space  charge  in  the  blocking  layer  and  at  the  active 
layer/blocking  layer  interface.  This  will  allow  for  simulation  of  and  interpretation 
of  C-V  data  outside  the  range  of  the  restricting  assumptions  described  above. 
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IV.  MODELING  RESULTS 


A.  LOW  TEMPERATURE  C-V  MODELING 

The  question  of  space  charge  in  the  blocking  layer  lends  itself  perfectly  to 
application  of  computer  simulation.  Whereas  in  actual  devices  space  charge 
cannot  be  directly  measured  on  a  point-by-point  basis,  the  model  calculates 
values  for  all  the  parameters  shown  in  Figure  6  for  finite  steps  through  the 
device.  It  is  therefore  a  simple  matter  of  extracting  and  plotting  the  data  to  see 
the  equilibrium  space  charge  distribution  in  a  given  device  configuration.  Figure 
9  shows  an  example  of  this  space  charge  profile  throughout  a  BIB  device.  The 
values  between  11  and  17  //m  are  plotted  in  blue  to  indicate  that  the  space 
charge  changes  sign  in  this  region.  Most  importantly,  the  data  is  continuous 
across  the  device. 
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Figure  9  Space  charge  profile  for  a  Si  BIB  with  a  0.05  V  Bias  (Table  1). 
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The  natural  point  to  begin  modeling  the  low  temperature  C-V  profiling 
method  is  therefore  to  solve  for  capacitance,  Equation  3.4,  using  the  full  spatially 
determined  space  charge  (Qsc)  for  a  given  bias.  This  means  we  no  longer  need 

to  tolerate  the  approximation,  Equation  3.3,  that  the  minority  doping  in  the  active 
layer  constitutes  all  the  space  charge.  Solving  for  capacitance  at  any  point  in  the 
device  requires  two  space  charge  profiles.  These  must  be  subtracted  to  solve  for 
a  change  in  space  charge  ( 8QSC )  over  a  small  change  in  bias  (dV).  In  the 

computer  simulations,  the  small  change  in  bias  was  kept  constant  at  0.1  mV. 
Subtracting  the  two  space  charge  profiles  gives  both  the  magnitude  and  the 
distribution  of  the  change  in  space  charge.  By  applying  this  method  for  multiple 
biases,  it  is  possible  to  probe  through  the  device  as  shown  in  Figure  10. 


1200 

1000 

800 

600 

400 

200 

0 

-200 


2 

a> 


o 

'Z 

-4-> 

o 

o 

LU 


Figure  10  Modulated  change  in  space  charge  (left  axis,  solid)  and  electric 
field  (right  axis,  dashed)  as  a  function  of  position  for  a  case  with 

an  ideally  sharp  interface. 
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Absorbing  Layer  Majority  Doping 

2.0  x  1017  cm'3 

Absorbing  Layer  Minority  Doping 

1 .0  x  1013  cm'3 

Absorbing  Layer  Thickness 

20  //m 

Blocking  Layer  Majority  Doping 

1 .0  x  1012  cm'3 

Blocking  Layer  Minority  Doping 

1.0  x  1011  cm'3 

Blocking  Layer  Thickness 

5  //m 

Temperature 

4  K 

Bohr  Radius 

1 5  x  1 0'8  m 

Contact  Grade  Parameter 

1 .0  x  10'5  cm 

Optical  Flux  Rate  y 

5.0  x  10'8  s'1  for  Photo-Current 

Table  1  Si  BIB  parameters  used  for  C-V  simulations 


The  idealized  case  in  Figure  10  has  the  device  parameters  shown  in  Table 
1,  with  a  sharp  grade  (g=2x10'7  cm)  between  the  active  layer  and  blocking  layer. 
Slight  variations  in  device  parameters  will  be  made  throughout  this  section,  but 
they  will  be  referenced  back  to  these  general  conditions.  In  Figure  10,  it  is 
apparent  that  the  change  in  space  charge  occurs  almost  entirely  at  the  outer 
boundary  of  the  depletion  region.  Integrating  this  space  charge  results  in  the 
total  change  in  space  charge,  and  this  can  be  inserted  into  Equation  3.4  to  solve 
for  the  capacitance.  For  the  range  of  bias  voltages  initially  selected  it  is  now 
possible  to  produce  a  plot  of  1/C2  versus  V  (Figure  11).  The  minority  doping 
concentration  is  proportional  to  the  inverse  of  the  slope  of  1/C2  as  seen  in 
Equation  3.6.  The  active  layer  minority  doping  level  shown  in  Table  1  for  this 
idealized  case  matches  the  calculated  values. 
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Applied  Bias  (Volt) 


Figure  11  1/C2  and  minority  doping  as  a  function  of  applied  bias. 

Parameters  listed  in  Table  1. 


One  adjustment  has  been  made  in  the  simulation  results  shown,  and  this 
is  due  to  capacitance  changes  in  the  immediate  contact  areas.  Free  carriers  in 
the  contacts  diffuse  into  the  device  and  recombine  over  a  given  distance  setting 
up  a  space  charge  region  known  as  a  Debye  tail.  The  space  charge  in  this 
region  also  changes  with  the  bias,  as  seen  in  Figure  12.  Due  to  its  higher 
concentration  of  carriers,  the  Debye  tail  will  have  a  lower  resistivity  than  the  rest 
of  the  device.  This  lower  resistivity  region  will  also  have  a  smaller  RC  time 
constant.  In  experiments  contact  effects  are  therefore  removed  with  other 
parasitic  effects  by  subtracting  a  higher  frequency  AC .  In  simulations,  the  model 
outputs  the  total  spatial  variation  of  the  change  in  space  charge,  so  this  high 
frequency  subtraction  can  be  accounted  for  by  setting  near  contact  values  equal 
to  zero. 
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Figure  12  Space  charge  variations  that  include  near  contact  space  charge 

(far  right).  Parameters  listed  in  Table  1. 

B.  INTRODUCING  INTERFACE  GRADES 

The  result  of  the  calculations  shown  in  Figure  1 1  demonstrates  that  low 
temperature  C-V  profiling  is  accurate  for  cases  that  match  the  assumption  of  a 
high  purity,  space  charge  free  blocking  layer.  The  focus  of  this  research  is  to 
investigate  the  effects  of  space  charge  in  the  blocking  layer  on  the  C-V  profiling 
method,  with  the  goal  of  determining  if  it  could  be  responsible  for  unexpected 
profiling  results.  Figure  13  shows  one  such  result  from  a  Raytheon  Vision 
Systems  profiling  run  on  material  for  Si:As  BIBs  [14]. 
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JWST  Assy  30748  Wafer  9406  Lot  1424-01  DieE09 
10.0  Hz  C-V,  Diode  D25, 7.66E-03cm2 


Figure  13  Low  temperature  C-V  profiling  results  (Courtesy  of  P.  Love). 

The  grade  between  doping  levels  at  the  active/blocking  layer  interface 
plays  a  significant  role  in  the  buildup  of  space  charge,  so  including  this  parameter 
in  the  modeling  is  very  important.  As  discussed  earlier,  one  of  the  difficulties  with 
growing  thin  pure  blockers  is  that  the  doping  levels  do  not  decrease  rapidly 
enough  from  heavily  doped  to  highly  pure.  Instead,  a  transition  region  exists 
where  the  material  properties  are  intermediate.  Figure  14  shows  one  such 
transition  region  as  measured  with  a  spreading  resistance  probe  for  a  Si:As  BIB. 
Starting  from  the  left  side  of  the  plot,  the  blocking  layer  has  a  doping  level  of 
2x1 013  cm'3  at  the  furthest  point  from  the  active/blocking  layer  interface,  but 
clearly  this  does  not  characterize  the  entire  blocking  layer.  In  fact,  2  //  m  into  the 
layer  the  doping  level  has  already  transitioned  to  IxlO14  cm'3.  This  characteristic 
must  be  taken  into  account  to  fully  understand  device  performance  [14]. 
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Figure  14  Spreading  resistance  measurement  showing  graded  region, 
blocking  layer  starting  from  the  left,  active  layer  -5-35  n  m 

(Courtesy  of  P.  Love). 

The  graded  region  is  modeled  in  the  program  as  previously  shown  in 
Figure  4,  allowing  variation  of  the  extent  of  this  transition  region.  The  first  step 
taken  which  included  the  grade  was  to  model  a  simplified  case  with  minimal 
space  charge  variation  between  the  two  regions.  In  this  case,  the  graded  region 
was  set  to  a  very  gradual  grade  parameter  of  g=2x10'5  cm,  and  the  majority 
doping  concentration  was  held  constant  at  2x1 017  cm'3  across  the  entire  device 
such  that  the  only  doping  variation  was  in  the  minority  doping  level.  Minority 
doping  levels  were  set  to  IxlO13  cm'3  in  the  active  layer  with  a  graded  transition 
to  5x1 012  cm'3  in  the  blocking  layer  with  all  other  parameters  the  same  as  shown 
in  Table  1.  Once  again  a  range  of  biases  were  chosen  and  a  set  of  two 
simulations  were  run  for  each  bias,  one  0.1  mV  higher  than  the  other.  Figure  15 
shows  that  the  compensation  calculations  closely  follow  the  expected  transition 
from  the  blocking  layer  value  of  5x1 012  cm'3  across  the  graded  interface  to  the 
active  layer  value  of  IxlO13  cm"3. 
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Applied  Bias  (Volt) 

Figure  15  Minority  doping  concentration  as  a  function  of  bias  for  g=2x10~5 
cm.  Inset  depicts  the  actual  position  of  space  charge,  and 
therefore  the  minority  doping  sample  area,  for  several  biases. 

C.  REALISTIC  DEVICE  PARAMETERS  WITH  GRADE 

The  previous  case  showed  that  the  low  temperature  C-V  profiling  method 
should  work  in  a  case  with  doping  variation  but  also  minimal  space  charge  in  the 
blocking  layer  and  at  the  interface.  We  now  repeat  the  simulations  from  the 
previous  section,  but  for  a  case  in  which  the  majority  doping  in  the  blocking  layer 
is  reduced  to  a  realistic  IxlO12  cm'3  and  the  minority  doping  is  reduced  to  IxlO11 
cm'3  (all  Table  1  values).  Changing  the  blocking  layer  majority  doping  alters  the 
space  charge  distribution,  and  the  role  of  the  grade  between  the  blocking  layer 
and  the  active  layer  becomes  significant.  Figure  16  and  Figure  17  demonstrate 
the  effects  of  varying  grade  parameters  on  C-V  results  for  otherwise  identical 
structures. 
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Figure  16  1/C2  as  a  function  of  applied  bias.  Parameters  listed  in  Table  1 

with  varied  interface  grade. 


Applied  Bias  (Volt) 

Figure  17  Compensation  calculations  from  C-V  profiling,  calculated  from 
Figure  16,  as  a  function  of  the  applied  bias.  Parameters  listed  in 
Table  1  for  varied  interface  grade. 
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For  the  case  with  the  sharpest  grade,  g=2x10'7  cm,  the  calculated  minority 
doping  concentration  as  seen  in  Figure  17  is  almost  exactly  the  actual  simulated 
doping  level  of  IxlO13  cm'3.  It  remains  consistent  at  this  level  as  the  depletion 
width  is  varied.  As  the  grade  is  increased,  the  calculated  values  from  the  initial 
part  of  the  curve  for  the  initial  bias  increase  steadily.  In  addition  the  plots  begin 
to  show  surprising  step  functions  and  variations  at  higher  biases.  These  results 
show  similar  characteristics  to  the  experimental  results  shown  in  Figure  13. 

It  is  clear  that  space  charge  in  the  blocking  layer  and  at  the  active 
layer/blocking  layer  interface  affects  calculations  of  minority  doping 
concentration.  Standard  interpretation  of  the  C-V  profiling  method  ignores  any 
space  charge  in  the  blocking  layer.  It  assumes  that  small  changes  in  the  applied 
bias  will  change  the  electric  field  profile  only  at  the  end  of  the  depletion  width  in 
the  active  layer.  This  is  a  safe  assumption  for  the  idealized  cases  initially 
discussed  which  include  uniform  layers,  layers  with  very  similar  doping  levels,  or 
layers  with  different  doping  levels  sharply  graded  at  the  interface.  However,  in 
non-ideal  cases,  particularly  when  a  grade  is  introduced,  this  is  no  longer  a  valid 
approximation  as  seen  in  the  strikingly  similar  results  of  the  simulations,  Figure 
17,  and  the  experimental  results,  shown  in  Figure  13.  The  next  step  is  to 
determine  why  the  space  charge  causes  such  non-linear  results. 

D.  ANALYSIS 

When  the  blocking  layer  ceases  to  be  ideal,  the  device  begins  to  behave 
like  two  separate  depletion  regions;  changes  in  the  electric  field  occur  both  at  the 
anticipated  point  near  the  end  of  the  depletion  width  in  the  active  layer  and  in  the 
blocking  layer  due  to  local  space  charge.  This  double  depletion  effect  is  most 
evident  as  the  grade  parameter  is  increased.  For  the  largest  grade  parameter, 
g=4x10'5  cm,  Figure  18  shows  the  space  charge  variation  in  both  layers.  Clearly, 
in  this  case  the  assumption  of  a  space  charge  free  blocking  layer  ceases  to  be 
valid.  The  leap  from  the  blocking  layer  space  charge  profiles  seen  in  Figure  18 
to  the  minority  doping  calculations  seen  in  Figure  17  is  complex  and  requires  an 
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in-depth  look  at  what  is  actually  happening  in  the  blocking  layer  and  at  the  active 
layer/blocking  layer  interface. 


Position  (micron) 

Figure  18  Change  in  space  charge  as  a  function  of  bias  for  grade  g=4x10~5 

cm. 


The  best  starting  point  is  an  analysis  of  the  change  in  space  charge  (8QSC) 
profiles  in  the  active  layer,  which  have  already  been  introduced  and  explored  in 
Figure  10.  Each  SQ^C  profile  is  a  result  of  the  corresponding  0.1  mV  change  in 

applied  bias  between  the  two  simulated  runs.  In  the  active  layer  the  total  SQsc 
between  each  set  of  biases  decreases  as  the  overall  bias  is  increased.  This 
occurs  because  the  relationship  is  nonlinear,  and  the  capacitance  ( dQsc/dV ) 

varies  as  a  function  of  l/Vv  (Equation  3.4).  The  0.1  mV  variation  no  longer 
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causes  as  large  an  increase  in  8QSC  at  higher  bias,  so  the  measured  capacitance 
decreases  with  bias,  and  1/C2  increases,  in  the  ideal  case,  linearly  with  V. 

This  is  crucial  to  understanding  the  calculation  of  the  minority  doping 
concentration.  It  requires  that  the  derivative  of  l/C2  be  positive  for  the  calculated 
minority  doping  to  be  positive  (Equation  3.6).  The  only  way  this  can  occur  is  if 
the  total  8QSC  decreases  as  the  applied  bias  increases.  One  ramification  of  this  is 
that  the  extracted  minority  doping  concentration  approaches  infinity  if  the  slope  of 
l/C2  transitions  from  positive  to  negative.  This  can  be  seen  in  Figure  16  and 
Figure  17  for  the  g=1x10'5  cm  grade.  The  slope  of  l/C2  approaches  zero  and 
becomes  negative  and  the  calculated  minority  doping  approaches  positive  infinity 
before  becoming  a  negative  (meaningless)  number. 

A  straightforward  physical  description  of  the  results  seen  in  Figure  16  and 
Figure  17  is  difficult  because  the  nonlinear  results  arise  from  multiple  interrelated 
factors.  These  factors  are 

a.  the  change  in  the  effective  bias  in  the  active  layer  resulting  from 
blocking  layer  space  charge, 

b.  the  change  in  total  space  charge  because  of  blocking  layer  space 
charge  magnitudes, 

c.  the  spatial  variation  of  the  graded  dopant  concentrations. 

These  factors  help  in  understanding  three  key  features  in  Figure  17: 

1.  variation  in  predicted  value  with  grade, 

2.  negative  slope  in  1/C2  leading  to  negative  calculated  doping  levels, 

3.  step  functions  in  calculated  doping  levels. 
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Figure  19  Grade  effects  on  initial  (0.05  V)  calculated  minority  doping. 

Starting  with  the  variation  in  the  predicted  value  as  a  function  of  grade  (1), 
Figure  19  shows  the  overall  increase  in  calculated  doping  level  as  the  grade 
increases,  for  the  lowest  bias  used.  Although  the  calculated  value  deviates 
further  from  the  actual  value  as  the  grade  is  increased,  the  effect  is  clearly 
nonlinear.  Figure  20  shows  the  expanded  field  profiles  of  several  of  these  runs 
at  the  interface  for  increasing  grade  between  the  active  and  blocking  layer.  The 
change  in  field  at  the  interface  is  a  function  of  the  increase  in  grade  between  the 
two  regions.  At  this  low  bias,  positive  space  charge  in  the  blocking  layer 
immediately  adjacent  to  the  active  layer  is  responsible  for  a  significant  field 
variation  at  the  interface. 

One  effect  of  this  positive  space  charge  is  to  increase  the  field  at  the 
interface,  thereby  increasing  the  overall  extent  of  the  depletion  width  into  the 
active  layer,  as  shown.  Adjusting  the  amount  of  space  charge  in  the  blocking 
layer  will  therefore  change  the  depletion  width,  resulting  in  an  overall  change  in 
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the  total  bias  dropped  in  the  active  layer.  The  slight  0.1  mV  increase  in  bias 
causes  a  change  in  space  charge  in  the  blocking  layer.  This  in  turn  changes  the 
overall  depletion  width  into  the  active  layer  from  what  it  would  have  been  for  a 
0.1  mV  increase  in  bias  with  no  blocking  layer  variation.  The  resultant  8QSC  in  the 

active  layer  behaves  as  if  a  different  bias  had  been  applied  (a).  This  changes  the 
measured  capacitance  from  what  it  would  have  been  without  the  blocking  layer 
space  charge,  and  explains  why  the  extracted  values  of  minority  doping 
concentration  at  low  biases  increase  for  increased  grades. 


Figure  20  Expanded  field  profiles  near  the  interface  at  0.05  V  bias,  varied 

grade  parameters. 

Returning  to  the  8Qsc  profiles  in  Figure  18  it  is  important  to  look  closely  at 

the  characteristics  in  the  blocking  layer,  detailed  in  Figure  21.  First,  there  is 
significant  negative  change  in  space  charge  for  the  lowest  bias  of  0.05  V 
(decrease  in  positive  space  charge).  This  negative  change  will  cause  an  overall 
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decrease  in  field  from  the  case  with  no  change.  At  the  next  bias,  there  is  still  a 
negative  change,  but  it  is  smaller  so  the  overall  decrease  in  field  will  be  smaller 
than  for  the  first  case.  The  result  is  that  the  change  in  space  charge  in  the  active 
layer  is  larger  than  it  would  have  been  without  the  blocking  layer  effects,  which 
causes  a  smaller  1/C2  slope  and  an  increased  minority  calculation.  These  effects 
are  nonlinear  with  grade  as  seen  in  Figure  19  and  explain  why  the  initial 
calculated  minority  concentrations  increase  as  the  grade  increases. 


Position  (micron) 

Figure  21  Change  in  space  charge  in  the  blocking  layer,  grade  g=4x10~5  cm. 

The  next  key  feature  in  Figure  17  is  the  rapid  increase  and  then  sign 
change  in  the  minority  doping  calculations  (2).  This  characteristic  occurs  for  the 
g=4x10'5  cm  case  in  Figure  17  and  correlates  to  the  transition  to  a  negative  1/C2 
slope  in  Figure  16.  As  previously  discussed,  a  negative  calculated  minority 
calculation  occurs  when  the  change  in  8QSC  becomes  positive,  so  it  is  expected 

that  8QSC  steadily  decrease  as  the  depletion  width  extends  further  into  the  device. 
In  Figure  17  the  g=4x10'5  cm  case  becomes  negative  at  0.8  V  as  a  result  of  this 


35 


second  factor,  the  absolute  magnitude  of  space  charge  in  the  blocking  layer  (b). 
Looking  at  the  space  charge  profiles  in  Figure  18  and  Figure  21,  it  can  be  seen 
that  the  d Qsc  in  the  blocking  layer  increases  with  increased  bias  once  it  becomes 

positive.  It  becomes  a  more  and  more  significant  portion  of  the  total  8QSC  in  the 

device  until  finally  at  0.8  V  it  represents  a  larger  increase  than  the  decrease  in 
the  active  layer,  and  the  net  change  is  positive.  This  represents  a  total 
breakdown  of  the  assumptions  made  in  the  standard  method. 

The  last  key  feature,  and  most  complex,  is  a  simulated  step  function  in  the 
calculated  minority  doping  (3).  This  feature  combines  the  two  factors  already 
discussed,  variation  of  the  effective  change  in  bias  in  the  active  region  (a)  and 
the  absolute  magnitude  of  the  space  charge  in  the  blocking  layer  (b),  with  the 
third  factor  (c).  This  last  factor  is  the  grade’s  significant  minority  dopant  variation 
(c),  which  can  also  be  seen  in  Figure  20.  It  changes  the  fields  from  straight  line 
approximations  as  in  the  sharpest  case  (g=2x10'7  cm)  to  the  wider  curves  pulled 
further  into  the  blocking  layer  for  the  most  gradual  case  (g=8x10'5  cm).  The  net 
result  is  another  level  of  complexity  and  nonlinearity  not  included  in  the  low 
temperature  C-V  profiling  approximation  of  minority  doping. 

The  majority  of  the  data  seen  in  Figure  17  represents  the  combined 
effects  of  all  of  the  factors  listed.  This  coupling  of  the  many  blocking  layer  space 
charge  effects  leads  to  the  nonlinear  results  seen,  and  variations  cannot  be 
easily  attributed  to  any  one  cause.  Most  importantly,  this  result  exposes  the 
dangers  of  ignoring  space  charge  in  the  blocking  layer  when  calculating  minority 
doping  concentrations,  particularly  when  graded  regions  are  introduced.  The 
combined  effects  of  blocking  layer  space  charge  on  the  depletion  width  and  on 
the  overall  value  of  8QSC  can  result  in  deviation  from  actual  values,  negative 

calculated  values,  and  step  functions  for  low  temperature  C-V  calculated  minority 
doping  concentrations  as  seen  in  Figure  17. 
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V.  COMPENSATION  IN  ALTERNATE  BIAS 


A.  ALTERNATE  BIAS  OPERATION 

Many  of  the  fabrication  issues  crucial  to  extending  wavelength  response 
into  the  mid  to  far-IR  are  tied  to  the  difficulty  of  growing  highly  pure,  extremely 
thin  blocking  layers.  When  these  conditions  are  not  simultaneously  achieved  in 
the  standard  configuration  (see  Chapter  I,  Section  C),  performance  is  affected. 
In  one  limit  an  extremely  thin,  but  impure  blocking  layer  would  not  provide 
sufficient  blocking  of  hopping  current  and  would  allow  a  high  dark  current.  In 
another  limit  an  ultimately  pure  blocking  layer  but  with  a  graded  transition  region 
from  heavily  doped  to  highly  pure  would  have  the  same  results  in  a  thin  layer 
limit.  Here  the  blocking  layer  needs  to  be  grown  thick  enough  to  isolate  the 
graded  region  from  the  active  layer  [11,12], 
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Figure  22  Electric  field  as  a  function  of  position  for  several  blocking  layer 
thicknesses.  Parameters  the  same  as  Table  1  except  active  layer 
minority  doping  5x1 012  cm'3,  temperature  4.2  K,  bias  0.75  V. 
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Growing  thicker  blocking  layers  is  not  a  desirable  solution  in  standard  bias 
because,  as  shown  in  Figure  22,  a  thicker  blocking  layer  results  in  significant  loss 
in  depletion  width  in  the  active  layer  and  loss  in  charge  collection  and 
responsivity  at  a  given  bias.  Keeping  the  blocking  layer  as  thin  as  possible 
reduces  both  this  loss  of  responsivity  and  the  range  of  space  charge  effects 
identified  through  the  earlier  C-V  analysis. 

One  proposed  solution  to  these  problems  is  to  operate  under  an  alternate 
bias  [15].  In  this  configuration  the  bias  is  applied  to  the  active  layer  contact 
rather  than  the  blocking  layer  contact.  With  an  n-type  BIB  this  requires  applying 
a  positive  bias  to  deplete  into  the  compensated  minority  sites  (NA)  as  seen  in 

Figure  23.  Here  the  photo-ionized  donor  sites  will  be  pushed  toward  the  blocking 
layer  interface  where  the  hopping  current  will  be  stopped,  and  the  electrons  in 
the  conduction  band  will  be  collected  at  the  active  layer  contact. 


Figure  23  Schematic  diagram  of  alternate  biased  n-type  BIB.  Red  line 

indicates  the  electric  field  distribution. 
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For  the  alternate  biased,  p-type  BIB  in  Figure  24  a  negative  bias  must  be 
applied  to  deplete  into  the  positive  compensated  donor  sites  (N+).  Here,  the 
electric  field  will  pull  ionized  holes  in  the  valence  band  toward  the  active  layer 
contact  and  push  NA  carriers,  the  hopping  current,  toward  the  blocking  layer 
where  they  will  be  stopped. 


IR 


Figure  24  Schematic  diagram  of  alternate  biased  p-type  BIB.  Red  line 

indicates  the  electric  field  distribution. 

This  alternate  bias  configuration  is  desirable  because  it  causes  the  electric 
field  to  deplete  through  the  active  layer  first  and  results  in  a  larger  depletion 
region  in  the  active  layer  as  seen  in  Figure  25.  This  effect  alone  is  desirable 
because  it  results  in  larger  photocurrents  in  ideal  blocking  layer  cases,  although 
dark  current  may  also  be  larger.  The  advantages  become  more  pronounced 
when  realistic  blocking  layer  parameters  for  longer  wavelength  BIBs  are  taken 
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into  account.  The  alternate  bias  allows  for  significantly  thicker  blocking  layers 
while  maintaining  the  depletion  region  in  the  active  layer. 


Position  (microns) 

Figure  25  Electric  field  as  a  function  of  position  for  standard  and  alternate 

biases.  Parameters  same  as  Table  1. 

The  primary  concern  in  this  configuration  is  ensuring  that  the  device  does 
not  saturate  with  dark  current  as  a  result  of  the  alternate  bias.  The  hopping 
current,  which  results  from  charge  flow  along  the  dopant  band,  is  most  critical 
and  can  be  calculated  in  the  model.  Figure  26  shows  the  ionized  donor  current 
as  a  function  of  position  for  both  standard  and  alternate  bias.  The  purity  of  the 
blocking  layer  prevents  the  completion  of  the  circuit  which  would  allow  the 
hopping  current  to  flow  continuously  through  the  dopant  band.  In  standard  bias 
(n-type  BIB)  it  is  stopped  because  the  highly  pure  blocking  layer  prevents  the 
injection  of  ionized  donors  N+D  from  the  blocking  layer  contact  into  the  active 

layer.  In  alternate  bias  the  same  is  true,  only  the  ionized  donors  are  being 
pushed  toward  the  blocking  layer  and  the  high  purity  of  the  blocking  layer,  and 
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the  associated  low  mobility  for  hopping  conduction,  prevents  these  carriers  from 
reaching  the  contact. 


Figure  26  Hopping  current  as  a  function  of  position  for  multiple  biases  and 
configurations.  Parameters  the  same  as  Table  1  except  active 
layer  minority  doping  5x1012  cm'3,  temperature  4.2  K. 


Also  of  concern  in  alternate  bias  is  the  diffusion  of  charge  carriers  from  the 
contacts.  In  standard  bias,  the  carriers  diffuse  from  both  contacts  and 
recombine,  setting  up  near-contact  space  charge  regions.  The  polarity  of  the 
applied  field  is  such  that  in  standard  bias  the  charge  carriers  are  pulled  from  the 
active  layer  contact  into  the  active  layer,  and  pushed  from  the  blocking  layer  into 
the  blocking  layer  contact  on  the  other  side  as  seen  in  Figure  27.  The  net  result 
is  that  the  significant  Debye  tail  occurs  in  the  heavily  doped  active  layer.  The 
extent  of  the  Debye  tail  is  determined  by  the  trapping  of  free  carriers  by  ionized 

sites  and  is  therefore  proportional  to  A/Nmin  In  alternate  bias,  the  injecting 
contact  for  free  electrons  is  switched  to  the  blocking  layer  as  a  result  of  the 
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switch  in  bias.  Because  the  minority  dopant  concentration  determines  the  rate  of 
recombination  and  thus  the  length  of  the  Debye  tail,  the  purity  of  the  blocking 
layer  requires  that  it  be  thick  enough  to  terminate  the  contact  diffusion  current 
prior  to  the  active  layer/blocking  layer  interface.  Otherwise,  the  hopping  current 
would  have  the  necessary  charge  source  and  path  to  complete  the  circuit,  raising 
the  dark  current  to  an  unacceptable  level. 
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Figure  27  Plot  of  free  carrier  concentration,  showing  Debye  tails  for 
standard  and  alternate  bias.  Parameters  the  same  as  Table  1 
except  active  layer  minority  doping  5x1 012  cm'3,  temperature  4.2 

K,  bias  0.5  V. 


B.  ALTERNATE  CONFIGURATION  IN  GERMANIUM 

The  longer  Debye  tail  and  associated  dark  current  in  alternate  bias  could 
be  a  serious  setback  in  the  standard  configuration.  However,  the  strength  of  the 
alternate  bias  lies  in  its  operational  independence  from  blocking  layer  thickness. 
The  potential  impact  of  employing  alternate  bias  for  longer  wavelength  response 
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in  Ge  and  GaAs  is  significant  because  it  relaxes  blocking  layer  thickness 
requirements.  For  the  parameters  shown  in  Table  2,  Figure  28  demonstrates 
that  devices  in  alternate  bias  have  almost  no  decrease  in  depletion  width  for  a 
blocking  layer  thickness  of  50  //  m.  When  these  thicknesses  are  modeled  in  the 
same  devices  in  standard  bias,  Figure  28  shows  the  results  are  very  small 
depletion  widths  in  the  active  layer  and  a  resulting  small  photocurrent. 


Absorbing  Layer  Majority  Doping 

1 .0  x  1016  cm'3 

Absorbing  Layer  Minority  Doping 

4.0  x  1012  cm'3 

Absorbing  Layer  Thickness 

30  //m 

Blocking  Layer  Majority  Doping 

2.0  x  1012  cm'3 

Blocking  Layer  Minority  Doping 

1.0  x  1011  cm'3 

Blocking  Layer  Thickness 

50  //m 

Temperature 

2.5  K 

Bohr  Radius 

37  x  10'8  m 

Interface  Grade  Parameter 

2.0  x  10'7  cm 

Contact  Grade  Parameter 

1 .0  x  10'5  cm 

Optical  Flux  Rate  y 

1.0x1  O'9  s'1  for  Photo-Current 

Table  2  Ge  BIE 

1  parameters 
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Figure  28  Electric  field  as  a  function  of  position  for  standard  and  alternate 

bias  in  Ge.  Parameters  listed  in  Table  2. 

However,  the  success  of  this  approach  is  highly  dependent  upon  the 
minority  doping  level  in  the  blocking  layer.  When  the  same  Ge  device  is  modeled 
in  alternate  bias  with  a  higher  minority  doping  level,  the  voltage  drop  occurs 
almost  exclusively  in  the  blocking  layer.  The  increase  in  resistivity  of  the  blocking 
layer  causes  an  increase  in  voltage  drop  that  becomes  significant  for  the  large 
thickness.  Figure  29  illustrates  this  field  distribution,  and  shows  that  devices 
require  minority  dopant  concentrations  at  around  or  below  IxIO^cm'3  to  deplete 
adequately  into  the  active  layer. 
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Figure  29  Electric  Field  as  a  function  of  doping,  Ge.  Parameters  listed  in 
Table  2,  varied  minority  doping,  alternate  biased  0.2  V. 


Once  again  minority  doping  concentrations,  in  this  case  in  the  blocking 
layer,  are  crucial  to  the  success  or  failure  of  any  attempts  to  build  the  devices. 
The  alternate  bias  configuration  offers  a  significant  relaxation  of  blocking  layer 
growth  requirements  with  performance  improvements,  assuming  these  doping 
levels  can  be  held  to  a  low  enough  level  in  a  reasonable  growth  thickness  (~50 
//m). 
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VI.  CONCLUSIONS 


A.  SUMMARY 

Current  attempts  to  extend  BIB  response  into  the  40-350  //m  range  using 
Ge  or  GaAs  have  been  slowed  by  fabrication  difficulties.  Despite  this,  the 
potential  of  these  materials  for  large  arrays  and  for  new  applications  in  the 
terahertz  regime  continue  to  drive  the  industry  to  arrive  at  solutions.  The 
application  of  alternate  bias  configuration  for  Ge  devices  shows  that  blocking 
layer  thicknesses  can  be  increased  significantly,  thereby  alleviating  one  of  the 
key  constraints  on  realizing  a  far-infrared  Ge  BIB.  The  minority  doping 
concentration  remains  crucial  and  the  alternate  bias  will  only  work  if  these  levels 
are  kept  low  enough.  In  this  work,  it  has  been  established  that  the  minority 
doping  levels  in  the  blocking  layer  must  be  kept  at  or  below  IxlO11  cm'3  for  Ge 
BIBs.  GaAs  devices  likewise  rely  upon  low  compensation  levels,  although  for 
this  material  the  difficulties  arise  in  the  active  layer  rather  than  the  blocking  layer. 
The  end  result  is  that  performance  of  BIBs  using  either  of  these  materials 
requires  control  of  compensation  levels  in  both  the  active  and  blocking  layers  of 
these  devices.  Problems  with  accurately  measuring  the  minority  doping 
concentration  exacerbate  the  already  difficult  situation. 

The  numerical  model’s  ability  to  solve  for  space  charge  throughout  BIB 
devices  made  it  ideal  for  analyzing  the  low  temperature  C-V  profiling  method  of 
determining  minority  doping.  Using  this,  we  have  demonstrated  the  important 
role  of  space  charge  in  the  blocking  layer  in  the  application  of  C-V  profiling  for 
measuring  minority  doping  concentration,  particularly  for  the  case  of  graded 
interfaces  in  BIB  devices.  Modeled  results  closely  resemble  experimentally 
achieved  results  and  explain  the  presence  of  step  functions  and  overall  variation 
from  expected  values. 
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B.  FUTURE  WORK 

This  thesis  provides  a  basic  explanation  of  unexpected  minority  doping 
variations  in  low  temperature  C-V  measurements.  The  next  step  needs  to  be 
determination  of  the  actual  grade  between  the  active  layer  and  the  blocking  layer 
for  the  experimental  results  presented  (Figure  13).  Spreading  resistance  data 
can  be  obtained,  and  this  will  test  our  assertion  that  the  unexpected  data  is  a 
result  of  the  grade  between  the  active  and  blocking  layers. 

This  thesis  has  highlighted  ramifications  of  the  limiting  assumption  of  a 
space  charge  free  blocking  layer,  however  the  question  of  correcting  this  problem 
has  not  been  addressed.  A  thorough  investigation  of  this  space  charge  needs  to 
be  conducted  to  isolate  exactly  what  effects  it  has  upon  calculations.  A  full 
understanding  of  these  characteristics  should  make  it  possible  to  modify  the 
assumptions  and  solve  for  the  actual  minority  doping  concentration. 

Finally,  fabrication  of  a  BIB  device  specifically  for  testing  and  validation  of 
the  alternate  bias  configuration  is  crucial  for  acceptance  by  industry.  The 
potential  for  far-infrared  response,  particularly  for  Ge  BIBs,  must  be  explored  in 
the  laboratory  with  thicker  blocking  layers  to  test  the  results  achieved  in 
simulations  at  the  Naval  Postgraduate  School. 
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